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Abstract 

We have calculated the 6s — 7s parity nonconserving El transition ampli- 
tude (Epnc) in cesium. This calculation has been performed with higher 
numerical accuracy than our 1989 calculation [V.A. Dzuba, V.V. Flambaum, 
and O.P. Sushkov, Phys. Lett. A 141, 147]. Also the Breit interaction 
has been included and the radiative corrections estimated. Our final re- 
sult is E PNC = 0.902^1 ± (0.7%]jiea B (Qw/N). This represents an im- 
provement in the accuracy of the calculation from the 1% error claimed 
in 1989. This result corresponds to a nuclear weak charge for Cs, Qw = 
-72.39(l ± 0.4%(exp) ± 0.7%(theory)) . We conclude that there is no signifi- 
cant deviation from the Standard Model value —73.09(3). 



PACS: 32.80.Ys, 11.30.Er, 12.15.Ji, 31.30.Jv 
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I. INTRODUCTION 



Measurements of parity nonconservation (PNC) in atoms provide a very useful test of 
the Standard Model. In order to extract the value of the nuclear weak charge Qw from ex- 
periment, atomic structure calculations must be performed. The most precise measurements 
and calculations have been performed for cesium. In 1989, calculations for cesium reached 
an accuracy of 1% 0,0] ■ At this time, the experimental accuracy was at a level of 2%. In 
1997, Wood et al measured the PNC amplitude in cesium to an accuracy of 0.35% [|[]. In 
light of new measurements of quantities relevant to the PNC amplitude (which are now in 
better agreement with the calculated values), Bennett and Wieman analyzed the accuracy 
of the calculations jl]|| and concluded that it is 0.4% Q. With this accuracy, the nuclear 
weak charge deviates from the Standard Model value by 2.5 a. 

The suggestion that the actual accuracy of the calculations is 0.4% immediately raised 
questions about values of small corrections which were neglected in It was shown by 

Derevianko H that the contribution of the Breit interaction to the PNC amplitude, Ep^c, 
is substantially larger than previous estimates and reduces the deviation from the Standard 
Model. Sushkov pointed out || that neglected radiative corrections to Ep^c must also be 
considered. 

In this work we have performed a complete calculation of the 6s — 7s PNC amplitude 
in cesium. We use the same method of calculation as that used in the 1989 work |J (this 
method is described in Section [TJ) but with much better numerical accuracy. The method 
which we call "many-body perturbation theory in the screened Coulomb interaction" was 
developed to treat the most important sequences of higher-order correlation diagrams in 
all orders. Note that this is an all-order technique which is not a version of the popular 
coupled-cluster technique. This method has proven to be very accurate for alkaline atoms 
and has produced very good results in a number of calculations. We believe that it can 
hardly be improved in terms of incorporating more correlation diagrams. However, since 
more computer power is available to us now compared to what we had in 1989 it is important 
to check the stability of the results with the improved numerical accuracy. In this work we 
have also calculated the Breit and radiative corrections to the PNC amplitude. We have 
also accounted for the small shift in Ep^c due to consideration of the neutron distribution. 

We list the results of our work below. 

Repeating the 1989 calculation with much higher numerical accuracy, we have obtained 
the same result: Ep^c = 0.908 (in units ieap^—Qw/N) x 10~ n ). From the analysis below, 
we have concluded that the accuracy of this value is about 0.5%. 

Our calculation of the Breit interaction corrects Ep^c by —0.0055. This is in agreement 
with §g§. 

Radiative corrections to the weak charge of order a have been calculated in || (see also 
||10||). However, there are important radiative corrections of order Za 2 and Z 2 a s ln 2 (A/i? n ) 
which have recently been calculated in [Tl[] (note that the latter correction is larger). Here 
A is the electron Compton wavelength, and R is the nuclear radius. Their contribution 
to Epnc is 0.004. This contribution originates from the radiative corrections to the weak 
matrix element due to the Uehling potential (recently this contribution was calculated in 



T2|j). We should add that there are other important contributions which have the same 



magnitude. In the sum-over-states approach (see Section p| ) they correspond to radiative 
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corrections to the energy intervals (-0.003) and El amplitudes (0.003). Parametrically, 
they are proportional to Z 2 a 3 ln(l/Z 2 a 2 ). There are also corrections to the weak matrix 
elements of order Z 2 a 3 \n(X/ R). Our estimate of these corrections shows that they have 
opposite sign to the Uehling potential contribution and the same magnitude. As a final 
result we suggest the following estimate for the contribution of the radiative corrections to 
E PNC : 0.000 ±0.004. 

In this work we use an improved charge distribution compared to that used in 1989. This 
shifts E PNC very slightly, by 0.0007. 

The neutron distribution in nuclei is slightly different from the well-studied charge dis- 
tribution. This produces a small correction to Ep^c- In this work we have found this 
correction to be —0.0018. This is in agreement with the result of Ref. |13 . 

This leaves us with Ep^c = 0.902 ± 0.7%. Combined with measurements of Ep^c/P 
from |3j and (3 = 27.00(6)a^ (the average value of the most accurate results from fl4Tj4fl 
and pill) this gives Q w = -72.39(l ± 0.4%(exp) ± 0.7%(theory)) which deviates by 

(l.0±0.4(exp)±0.7(theory))% from the Standard Model value Q w = -73.09(3) p|. From 
our point of view, this does not look like a significant deviation from the Standard Model. 



II. METHOD OF CALCULATION 



The method that we use here was developed in the works |l|, ^7|JT8|| . As was emphasized 



in the introduction, it is an all-order technique, in terms of treating correlations, and is not 
a version of the popular coupled-cluster method. The dominating sequences of higher-order 
correlation diagrams correspond to real physical phenomena like screening of the Coulomb 
interaction and the hole-particle interaction. They are included in all orders in our tech- 
nique. While many-body perturbation theory in the residual Coulomb interaction does not 
converge, there is good convergence of our method "many-body perturbation theory in the 
screened Coulomb interaction" . Another strong point of the method is that its complexity 
does not go beyond the calculation of energies. The most complicated and time consuming 
part of the method is the calculation of the correlation potential S. £ is used to calculate 
single-electron Brueckner orbitals. The calculation of matrix elements with Brueckner or- 
bitals is easy and they already include most of the correlations. This makes the calculation 
of hyperfine structure, PNC, etc. as simple as the calculation of energies. 

The nuclear spin-independent weak interaction of an electron with the nucleus is 

G F 

H w = ^-^p(r)Q w j 5 (1) 

where Gp is the Fermi constant, Qw is the weak charge of the nucleus, 75 is a Dirac matrix, 
and p(r) is the nuclear density (note that because the nuclear weak charge in the Standard 
Model is approximately equal to the number of neutrons, the density p(r) should be taken 



as the neutron density; see Section pi]) . 

There are essentially two different methods by which the PNC El amplitude can be 
calculated: from a "mixed-states" approach (in which the external fields are taken into 
account in the electron orbitals) and from a "sum-over-states" approach (perturbation theory 
sum over intermediate opposite parity states). The most complete calculation using the 



3 



mixed-states method was performed in the work |IJ. In this work all-orders summation of 
dominating diagrams in the residual electron-electron interaction was included. This is the 
method used in the current work. (For the sum-over-states method, see the discussion in 
Section 0.) 

The method of calculation we use is applicable for iV-electron atoms with one valence 
electron. This method is particularly effective for cesium, compared to other heavy atoms, 
as in this case the external electron has very little overlap with the tightly-bound core, 
enabling the use of perturbation theory in the calculation of the residual interaction of the 
external electron with the core. 

The calculations start from the relativistic Hartree-Fock (RHF) method in the V^ -1 
approximation. The single-electron RHF Hamiltonian is 

H = co • p + (P - l)c 2 - Za/r + V^ 1 , (2) 

a and /3 are Dirac matrices and p is the electron momentum. For cesium, the accuracy of 
the RHF energies is of the order of 10%. 

We use the time-dependent Hartree-Fock (TDHF) method (which is equivalent to the 
random-phase approximation with exchange) to calculate the interaction of external fields 
with atomic electrons. In this paper we deal with two external fields: the electric field 
of the photon (El transition amplitudes) and the weak field of the nucleus. In the RHF 
approximation the interaction between an external field H ext and atomic electrons is defined 
by the matrix element (ip2\H ext \4 ! i) y where ip\ and ip2 are RHF orbitals. Inclusion of the 
polarization of the atomic core by an external field is reduced to the addition of a correction 
SV (which is the correction to the Hartree-Fock potential due to the interaction between the 
core and the external field) to the operator which describes the interaction, (V^l-^ext+^l^i)- 

The TDHF contribution to Ep^c between states 6s and 7s in the mixed-states approach 
is given by 

E T P ^ F = (^ 7 s\H E1 + 5V E1 \5^ S ) + (i> 7s \H w + 5V w \X 6s ) + ^ 7s \5V Elw \^ s ) (3) 

(of course, this term can instead be expressed in terms of corrections to if>? s )- Here 5ip and 
SVw denote corrections to single-electron RHF wavefunctions and the Hartree-Fock potential 
caused by the weak interaction. These corrections are found by solving 

(H -e)61> = -(H w + 5fr w )il> . (4) 

The positive (negative) frequency corrections X (Y) due to the El field of the photon are 
found from the equations 

(H -e- u)X = -{Hex + SV E1 )if> (5) 
(H -e + u)Y = -(H E1 + 5V E1 )^ . (6) 

The correction to the potential due to the El field is 5V E i (SV E1 ). The correction SV E iw to 
the core potential is due to the simultaneous action of the weak field and the electric field 
of the photon. 

The TDHF contribution to E PNC (Eq. |3|) corresponds to calculating the lowest-order 
PNC diagrams presented in Fig. [I] with the core polarization diagrams (see Fig. 0) included 
in all orders in the Coulomb interaction. 
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The wavefunctions are improved by taking into account electron-electron correlations. 



The correlations are calculated using the "correlation potential" method [|19[ which corre- 
sponds to adding a non-local correlation potential E to the potential V^" 1 in the RHF 
Hamiltonian @ and then solving for the states of the external electron. The correlation 
potential is defined such that its average value coincides with the correlation correction to 
energy, 5E{ = (ifti\Y}\iJ)i) . The correlation potential is calculated by means of many-body 
perturbation theory in the residual Coulomb interaction 

N N i N 

i=l i<j I * 3 1 i=l 

where H is the exact Hamiltonian of an atom. The lowest-order correlation diagrams 
(second-order in U) are presented in Fig. [3[ 

Brueckner orbitals are obtained by adding S^ 2 ^ to the Hartree-Fock potential V N ~ 1 and 
solving the equation 

(H + S (2) - e)V = (8) 

iteratively for the states of the external electron. This corresponds to chaining the self- 
energy operator to all orders (Fig. This chain of diagrams is enhanced by the small 
denominator, corresponding to the excitation energy of an external electron (in comparison 
to the excitation of a core electron). At this level of calculation (with "bare" E) the accuracy 
for energy levels for Cs is about 1%. 

Using the correlation potential method and the Feynman diagram technique we include 
into two series of dominating higher order diagrams which are calculated in all orders 



of perturbation theory [pj, p!7|,[18| . These are screening of the electron-electron interaction 
(Fig. |5|) and the hole-particle interaction (Fig. |6|). The electron-electron screening is a 
collective phenomenon. The corresponding chain of diagrams is enhanced by a factor ap- 
proximately equal to the number of electrons in the external closed subshell (the 5p electrons 
in Cs). The hole-particle interaction is enhanced by the large zero multipolarity diagonal 
matrix elements of the Coulomb interaction. This interaction accounts for the discrete spec- 
tra in noble gases. We will denote the dressed self-energy operator by E (Fig. 0). With this 
E the Brueckner energies for Cs have an accuracy of the order of 0.1%. 

The wavefunctions can be further modified by placing a coefficient before E such that 
the corresponding energy coincides with the experimental value. This fitting of the Brueck- 
ner orbitals can be considered as a way of including other higher-order diagrams into the 
calculations. 

If we use Brueckner orbitals instead of RHF orbitals to calculate the PNC amplitude 
in Eq. ^, then we include all-orders in E contributions to Ep^c- However, the correlation 
potential is energy- dependent, E = S(e). So, in first order, we should consider the proper 
energy dependence. The first-order in E corrections to Ep^c are presented diagrammatically 
in Fig. [S[ We can write these as 

(i> 7s \£(e 7s )\5X 6s ) + (8ip 7s \£(e 7s )\X 6s ) + (5Y 7s \t(e 6s )\^ s ) + (Y 7s \£(e 6s )\5^ s ) . (9) 

The non-linear in E contribution can be found by subtracting from the all-orders result the 
first-order value found in the same method. 
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The correlation corrections to Ep^c we have considered so far are usually called 
"Brueckner-type" corrections. (In this case the external field interacts with the external 
electron lines.) There are also contributions to Ep^c i n which the external field acts inside 
the correlation potential (see Fig. ||). Those diagrams in which the El interaction occurs 
in the internal lines are known as "structural radiation", while those in which the weak 
interaction occurs in the internal lines are known as the "weak correlation potential" . There 
is another second-order correction to the amplitudes which arises from the normalization of 



states fL9|| . The structural radiation, weak correlation potential, and normalization contri- 
butions are suppressed by the small parameter E ext /E core ~ 1/10, where E ext and E int are 
excitation energies of the external and core electrons, respectively. 

In this work we also calculate the contribution of the Breit interaction. We do not discuss 
the method of calculation here, but refer the interested reader to the works |j?|,|20" . 



We postpone the discussion of our method for the calculation of radiative corrections to 
Section [[Vj. 



III. 65 - 7S PNC AMPLITUDE 

The results of our calculation for the 6s — 7s PNC amplitude are presented in Table [[[ No- 
tice that the time-dependent Hartree-Fock value gives a contribution to the total amplitude 
of about 98%. The point is that there is a strong cancellation of the correlation corrections 
to the PNC amplitude. The stability of the PNC amplitude compared to other quantities 
in which the correlation corrections are large will be discussed in more detail in Section |V]. 
Notice that the values are in agreement with our 1989 result 0.908 x 10 _11 iea^(— Qyy/N) 
(see "Subtotal" of Table | for the current calculation). The higher numerical accuracy of 
the current work has therefore not changed the previous result. 

The mixed-states approach has also been performed in and || to determine the 
PNC amplitude in cesium. However, in these works the screening of the electron-electron 
interaction was included in a simplified way. In [Q] empirical screening factors were placed 
before the second-order correlation corrections S^ 2 ^ to fit the experimental values of energies. 
Kozlov and Porsev introduced screening factors based on average screening factors calculated 
for the Coulomb integrals between valence electron states. The results obtained by these 
groups (without the Breit interaction, i.e., corresponding to the Subtotal of Table [I]), 0.904 |2| 
and 0.905 ||, are close to ours. To be sure that we understand this difference, we performed 
a pure second-order (i.e., using £( 2 )) calculation and fitted the energies (as was done in [0) 
and reproduced their result, 0.905. 

In the work a calculation using the sum-over-states method was also performed. In 
the sum-over-states approach the 6s — 7s PNC amplitude is expressed in the form 



_ y> A7s\H w \np)(np\H E1 \6s) (7s\H E1 \np)(np\H w \Qs) \ 
pnc — 2^ V E- — e Ea - E ' 



J ls ■ 1 - J np - £ - / 6s J-'np 



The authors of reference include single, double, and selected triple excitations into their 
wavef unctions. Note, however, that even if wavefunctions of 6s, 7s, and intermediate np 
states are calculated exactly (i.e., with all configuration mixing included) there are still some 
missed contributions in this approach. Consider, e.g., the intermediate state 6p = 5p 6 6p. It 
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contains an admixture of states 5p 5 ns6d: 6p = 5p 6 6p + a5p 5 ns6d + .... This mixed state is 
included into the sum fllOD. However, the sum (|I~0"D must include all many-body states of 



opposite parity. This means that the state 5p 5 ns6d = 5p 5 ns6d — a5p^6p + ... should also 
be included into the sum. Such contributions to E PNC have never been estimated directly 
within the sum-over-states approach. However, they are included into the mixed-states 
calculation. The result of the sum-over-states approach, 0.909, is very close to the result 
of the mixed-states aproach, 0.908. It is important to note that the omitted higher-order 
many-body corrections are different in these two methods. This may be considered as an 
argument that the omitted many-body corrections in both calculations are small. Of course, 
here we assume that the omitted many-body corrections to both values (which, in principle, 
are completely different) do not "conspire" to give exactly the same magnitude. 

Therefore we will take 0.908 for the value of E PNC (Subtotal of Table |) as this corre- 
sponds to the most complete mixed-states calculation and is in agreement with the sum- 
over-states calculation of reference 0. 

With Breit, our result becomes 0.903 x 10 _11 iea^(— Qw/N). This correction is in agree- 
ment with 

We use the two-parameter Fermi model for the proton and neutron distributions: 

p(r) = p [l + exp[(r - c)/o]] , (11) 

where t = a(41n3) is the skin-thickness, c is the half-density radius, and po is found from 
the normalization condition / pdV = 1. In 1989 the thickness and half-density radius for 
the proton distribution were taken to be t p = 2.5 fm and c p = 5.6149 fm (corresponding to 
a root-mean-square (rms) radius (r 2 ) 1 / 2 = 4.836 fm). In this work we have used improved 
parameters t p = 2.3 fm and c p = 5.6710 fm ((t" 2 ) 1//2 = 4.804 fm) pl| . This changes the 
wavefunctions slightly, leading to a very small correction to the PNC amplitude of 0.08% 
(0.0007). (This is in agreement with a simple analytical estimate: the factor accounting for 
the change in the electron density is ~ (4.804/4.836)- z a ~ 0.1% .) In the work § we used 
the proton distribution in the weak interaction Hamiltonian (Eq. Q). In the current work we 
have found the small correction to Ep^c which arises from taking the (poorly understood) 
neutron density in Eq. [I]. We use the result of Ref. [22] for the difference Ar np = 0.13(4) fm 



in the root-mean-square radii of the neutrons (r 2 ) 1//2 and protons (?"p) 1//2 . We have considered 
three cases which correspond to the same value of (r 2 ): (i) c n = c p , a n > a p ; (ii) c n > c p , 
a n > a p ; and (iii) c n > c p , a n = a p (using the relation (r 2 ) w |c 2 + |7r 2 a 2 ). We have found 
that Epmc shifts from —0.18% to —0.21% when moving from the extreme c n = c p to the 
extreme a n = a p . Therefore, Ep^c changes by about —0.2% (—0.0018) due to consideration 
of the neutron distribution. This is in agreement with Derevianko's estimate, —0.19(8)% 

In the next section we discuss the radiative corrections to Ep^c- Our estimate of these 
corrections does not shift the PNC amplitude. 
Therefore, we have 

E PNC = 0.902 x l0~ n iea B (-Q w /N) (12) 

as our central point for the PNC amplitude. The error will be estimated in the following 
sections. 
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IV. QED-TYPE RADIATIVE CORRECTIONS TO ENERGY LEVELS, 
WAVEFUNCTIONS, AND THE PNC AMPLITUDE 



The radiative corrections to the weak charge Qw have been calculated for the free elec- 
tron. However, an electron in a heavy atom is bound, and this produces additional radiative 
corrections proportional to a(Za) n , n = 1,2, .... Recently such corrections were considered 
by Milstein and Sushkov fill . They found that the most important are corrections enhanced 
by the large parameter ln(A/i2), where A = h/mc is the electron Compton wavelength and 
R is the nuclear radius. This type of correction arises from the radiative corrections to 
the electron wavefunction near the nucleus. In this region the s-wave and pi/2-wave (lower 
Dirac component) electron densities are singular, |?/>(r)| 2 ~ r -z' 2 a 2 _ rp^g rac ii a ti ve correc- 
tions modify the potential at small distances r < A, V(r) = —Za(l + 5)/r. Correspondingly, 
the electron wavefunctions change, |?/>(r)| 2 ~ r~ z a ( 1+<5 ) for r < A. This gives the radiative 
correction factor for the electron density inside the nucleus, 

Hf~(^) ZVM = -p(^Vln( Vfl )). (13) 
For the Uehling (vacuum polarization) potential 5 ~ a ln(A/r) |[23|| . This gives an additional 



power of the large parameter \n(X/R). This leads Milstein and Sushkov UTTJ to conclude that 
the Uehling potential gives a dominating radiative correction to Epnc, ~ Z 2 a 3 \n 2 (X/R). 
Numerical calculations of the Uehling potential contribution have been performed in fl2fl 
and in the present work. This radiative correction increases Ep^c by 0.4%. 

Milstein and Sushkov [11] demonstrated that there are no other radiative corrections 
which are enhanced by ln 2 (A/i?). However, any correction to the potential with nonzero 
5(R) ~ a gives a correction to the electron density ~ Z 2 a 3 ln(X/R). We demonstrate below 
that such corrections are also important and give a contribution of opposite sign to that of 
the Uehling potential. 

Let us start our discussion from the radiative corrections to energy levels (the Lamb 
shift). The calculation of the shift can be divided into two parts: one in which the electron 
interaction with virtual photons of high-frequency are considered, and one in which virtual 
photons of low-frequency are considered. 

In the high-frequency case the external field (the strong nuclear Coulomb field) need 
only be included to first order. In this case the contributions to the Lamb shift arise from 
the diagrams presented in Fig. [10]. The contribution of the Uehling potential (Fig. |10"l(a)) 
to the Lamb shift is very small. The main contribution comes from the vertex correction 
(Fig. 1101(b)). (In the case of a free electron the vertex diagrams give the electric f{q 2 ) 



and magnetic g(q 2 ) formfactors.) The perturbation theory expression for f(q 2 ) contains 
an infra-red divergence and requires a low- frequency cut-off parameter k - see, e.g., p3| . 



Assuming q 2 -C m 2 c 2 , the high-frequency contribution to the Lamb shift can be presented 
potential given by the following expression [B3] 



5$(r) 



6$f + 5$ 



u 



5% 



a ^ 2 /, m 11 1\ , . ah , lS 



3nm 2 c 2 ^ 2k 24 h' 1 ' Airmc 
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For the Coulomb potential, A$ = —4iTZe5(r). Here the last long-range term (S$> g ) comes 
from the anomalous electron magnetic moment (g(0)); the infra-red cut-off parameter k 
appears from the electric formfactor f{q 2 )- This term with the large In ^ gives the dominant 
contribution to the Lamb shift of s-levels. The infra-red divergence for k — > indicates the 
importance of the low-frequency contribution for this term. 

If we go beyond the approximation q 2 « m 2 c 2 , the term 5$/ should be associated 
with a non-local self-energy operator £ ra d(r, r', E) with typical values |r — r'| < and 
t ~ r' < — . We need this operator in a simple limit, E « mc 2 . In this case we can 
approximate it by a two-parametric (A, b) potential 

5 $ = -A -e~ 6 ^$ . (15) 

7T 

The parameters A and b in 5$/ can be found from the fit of the Lamb-shift of the high 
Coulomb levels 3s, 4s, 5s and 3p, 4p and 5p (in one-electron ions) which were calculated 
as a function of the nuclear charge Z in Refs. |[24|| . We have checked that A = 1.25 and 
b = 1 fit all these Lamb shifts quite accurately (we have found slightly different potential 
strengths for s- and p-waves, A s = 1.17 and A p = 1.33). The anomalous magnetic moment 
contribution is 

5$ -_j_^?_ 7 .v$. (16) 
Anmc 



The Uehling potential for a finite nucleus is given by p5J (in atomic units (h = m = e = 1, 

a = 1/c)) 

f ** = ~ jf & *,(*) f dt v^T(i + ^) (e~ 2 ^ - e- 2 ^) , (17) 

where p(x) is the nuclear charge density. It is more convenient to use a simpler formula for 
5$u f° r r > R, R is the nuclear radius, 



*Mr) = m^jpfdt v^I(i + ^)e^/(a;) , (18) 
I( x ) = - e * + e ~ x + xe x + xe~ x , x = 2tR/a , (19) 

and take 5&u{ r < R) — $&u{ r — R)- There is practically no loss of numerical accuracy in 
this approximation since a typical scale for the variation of 5&u(r) is given by the electron 
compton length — >> R. The radiative corrections for Cs energy levels are presented in 
Table 0. 

Note that for the most important term 5$/ we do not use the assumption Za << 1 
since we fitted the exact results for the single-electron ions. The contribution of the Uehling 
potential 5§u to the Lamb shift is always small. Also, Milstein and Strakhovenko have shown 
in Ref. |26j that higher Za corrections are numerically not important for this potential. The 



potential <5$ 9 due to the magnetic formfactor is a long-range one. This also hints that there 
are no large higher Za corrections here. 

We can use 5$ to estimate the contribution of QED-type radiative corrections to the 
electron wavefunction and PNC amplitude Ep^c- Note that it is not enough to calculate 
the radiative corrections to the matrix element of the weak interaction (n'pi/2\Hw\ns) . 
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Corrections to the energy intervals like 6s — Qp are also important since these intervals are 
small at the scale of the atomic unit (~ 1/20) and sensitive to perturbations. The change in 
the energies also influences the large- distance behavior of the electron wavefunctions which 
determine the usual El amplitudes in the sum-over-states approach. Thus, the simplest 
way to proceed is to include 5$ into the Dirac-Hartree-Fock equations and then perform all 
calculations. 

The results are the following: the Uehling potential 5$u increases Ep^c by 0.41% (in 
agreement with ]12[); the contribution of <5$ g is small, —0.03% (due to cancellation of the 
contributions of the corrections to the s-wave and p-wave); and the contribution of 5&f is 
—0.65%. The total result of 5$ is —0.27% (and opposite to that of the Uehling potential). 

Note that the replacement of £ ra d(r, r', E = 0) by a local parametric potential (5$/) 
may be quite a crude approximation. Therefore, we have tested the sensitivity of the result 
to variation of the parameter b (parameter A is a function of b, it is found from the fit of the 
Coulomb energy levels). The increase of the radius of the potential 5Qf two times (b = 0.5) 
increases E PNC by 0.3% (the total contribution of 5$ is 0.04%). 

We do not discuss here the Za correction to the Z-boson exchange vertex. Milstein and 
Sushkov [JTTJ] have shown that this correction does not contain the large parameter ln(A/i?) if 
the calculations are performed in the Landau gauge and they concluded that this correction 
should be small. 

Here we should note that the radiative corrections to the wavefunctions (as well as the 
wavefunctions themselves) are not gauge invariant. However, the change of the electron 
density inside the nucleus is a gauge invariant, observable phenomenon. This makes the 
consideration presented above meaningful. 

Of course, our estimate of QED-type radiative corrections to E PNC should not be consid- 
ered as an accurate calculation. The aim is to show that the consideration of just the single 
Uehling potential contribution to the radiative corrections can give the wrong impression. 
Indeed, this contribution to Ep^c is enhanced by the large parameter ln 2 (A/i?). However, 
this only makes it comparable to other contributions which in the case of energy levels were 
an order of magnitude larger than that of the Uehling potential. A conservative estimate of 
the radiative correction contribution can be presented as 0.0 ± 0.4%. This range covers our 
two values (-0.27% and 0.04%) as well as the value 0.4% obtained in [0,[12|. 



V. ESTIMATE OF ACCURACY OF PNC AMPLITUDE 

We have estimated the error of the PNC amplitude in a number of different ways. There 
are two main methods: (i) root-mean-square (rms) deviation of calculated energy inter- 
vals, El amplitudes, and hyperfine structure (hfs) constants from the accurate experimental 
values; (ii) influence of fitting of energies and hyperfine structure constants on the PNC 
amplitude. 



A. Root-mean-square deviation 

Remember that the PNC amplitude can be expressed as a sum over intermediate states 
(see formula [L0|). Each term in the sum is a product of El transition amplitudes, weak 
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matrix elements, and energy denominators. There are three dominating contributions to 
the 6s - 7s PNC amplitude in Cs [§: 

(7s\H E1 \6p)(6p\H w \6s) (7s\H w \6p) (6p\H E1 \6s) (7s\H E1 \7p) (7p\H w \6s) 
hpNC ~ p p + p p + p p + ••• 

= -1.908 + 1.493 + 1.352 + ... = 0.937 + ... . (20) 

While we do not use the sum-over-states approach in our calculation of the PNC amplitude, it 
is instructive to analyze the accuracy of the El transition amplitudes, weak matrix elements, 
and energy intervals which contribute to Eq. |2(] as they have been calculated using the same 
method as that used to calculate E PNC - 

Let us begin with the energy intervals. The calculated ionization energies are presented in 



Table |TJ. The Hartree-Fock values deviate from experiment by 10%. Including the second- 
order correlation corrections reduces the error to ~ 1%. When screening and the hole- 
particle interaction are included into E( 2 ) in all orders, the energies improve, ~ 0.1%. The 
percentage deviations from experiment of the energy intervals of interest are: Eq s —Eq p , —0.3; 
E 7s — Eq p , 0.8; and Eq s — E 7p , —0.01. The rms error is 0.5%. We can in fact reproduce energy 
intervals exactly by placing coefficients before the correlation potential. Because this fitting 
of the energies appears to improve the wavefunctions (e.g., electromagnetic amplitudes and 
hyperfine structure constants improve) we will use this procedure in the following estimates. 

The relevant radial integrals (El transition amplitudes) are presented in Table [IV]. These 
were calculated with the energy-fitted "bare" correlation potential X^ 2 - 1 and the (unfitted and 
fitted) "dressed" potential E. Structural radiation and normalization contributions were also 
included. In Table [Vj the percentage deviations of the calculated values from experiment 
are listed. Without energy fitting, the rms error is 0.3%. Fitting the energy improves the 
accuracy, ~ 0.1%. 

We cannot directly compare weak matrix elements with experiment. Like the weak 
matrix elements, hyperfine structure is determined by the electron wavefunctions in the 
vicinity of the nucleus, and this is known very accurately. The hyperfine structure constants 
calculated in different approximations are presented in Table |VT[ Corrections due to the 
Breit interaction, structural radiation, and normalization are included. The percentage 



deviations from experiment are shown in Table |V11|. The rms deviation of the calculated 



hfs values from experiment using unfitted E is 1.1%. With fitting, the rms error in the 
pure second-order approximation is 0.3%; with higher orders we get 0.6%. We are, however, 
trying to estimate the accuracy of the s — p weak matrix elements. It makes more sense 



for us to use the square-root formula, ynfs(s)hfs(p). The errors are presented in Table [VTT . 
Notice that by using this approach the error is much smaller. Without energy fitting, the 
rms error is 0.4%. With fitting, the rms error in the second-order calculation (S^) and full 
calculation (E) is 0.2%. 

From this section we can conclude that the rms error for the relevant parameters is 
somewhere in the range 0.2 — 0.5%. 

Note that from this analysis the error for the sum-over-states calculation of Ep^c would 
be larger than this, as the errors for the energies, hfs constants, and El amplitudes contribute 
to each of the three terms in Eq. [H]. However, in the mixed-states approach, the errors do 
not add in this way. We get a better indication of the error of our calculation of Ep^c in 
the next section. 
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B. Influence of fitting on the PNC amplitude 



In the section above we presented calculations in three different approximations: with 
unfitted £, and with T,^ and £ fitted with coefficients to reproduce experimental ionization 
energies. The errors in these approximations are of different magnitudes and signs. We now 
calculate the PNC amplitude using these three approximations. The spread of the results 
can be used to estimate the error. 

The results are listed in Table [EXI. It can be seen that the PNC amplitude is very stable. 



The PNC amplitude is much more stable than hyperfine structure. This can be explained 
by the much smaller correlation corrections to PNC (compare Table | with Table |V^) . The 
stability of E PNC may be compared to the stability of the usual electromagnetic amplitudes 
where the error is very small (even without fitting). 

We have also considered the fitting of hyperfine structure using different coefficients 
before each £. Using the resulting wavefunctions, the PNC amplitude increased by 0.5%. 
This is the maximum deviation we have obtained. We will therefore use this as the estimate 
for the accuracy of the E PNC calculation. 

VI. CONCLUSION 

We have obtained the result 

Ep NC = 0.902(6) x W- n iea B (-Q w /N) (21) 

for our calculation of the 6s — 7s PNC amplitude in Cs. This is in agreement with other PNC 
calculations, however we would like to emphasize that our calculation is the most complete. 
The most precise measurement of the 6s — 7s PNC amplitude in Cs is 

- lm{E " Nc) = 1.5939(56)^ , (22) 
p cm 

where (3 is the vector transition polarizability. For (3 we use the value 

(3 = 27.00(6)4 ( 23 ) 



which is the average value of the most accurate results ^MM and [njJ3|. Using the conversion 
\e\/a 2 B = 5.1422 x 10 12 mV/cm, we therefore obtain for the weak charge of the Cs nucleus: 

Q w = -72.39(29) cxp (51) theory , (24) 

where the experimental error is obtained by adding in quadrature the error for (3 and the 
error for Im(EpNc)/ (3- This result deviates by (l.O ± 0.4(exp) ± 0.7(theory)^% from the 
Standard Model value Qw = —73.09(3) |16|. This does not represent a significant deviation. 
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TABLES 



TABLE I. Contributions to the 6s — 7s Ep^c amplitude for Cs in units 10 ll ieaB{—Qw/N). 
(S corresponds to the (unfitted) "dressed" self-energy operator.) 



TDHF 0.8885 
(^ 7s \t\5X 6s ) 0.0705 
{Sip 7s \t\X 6s ) 0.1857 

{5Y 7s \t\^ 6s ) -0.0760 

(Y 7s \£\5ip 6s ) -0.1420 

Nonlinear in £ correction -0.0198 
Weak correlation potential 0.0038 
Structural radiation 0.0025 

Normalization -0.0049 

Subtotal 0.9084 

Breit -0.0055 
New proton distribution 0.0007 

Neutron distribution -0.0018 
Radiative corrections 0.0 

Total 0.902 



TABLE II. Radiative corrections to RHF ionization energies; units —cm 1 . (See also Table PH.) 

6s 7s 6pi/ 2 7pi/ 2 

-18.4 -5.0 0.88 0.31 



TABLE III. Ionization energies for Cs in units —cm . 



State 



RHF 



S(2) 



Experiment 



6 s 
7s 

6P1/2 
7pi/ 2 



27954 
12112 
18790 
9223 



32415 
13070 
20539 
9731 



31420 
12863 
20276 
9657 



31407 
12871 
20228 
9641 



i Taken from [27| 
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TABLE IV. Radial integrals of El transition amplitudes for Cs in different approximations. 
The experimental values are listed in the last column, (a.u.) 



Transition 


RHF 


TDHF £W 

with fitting 


S 


t 

with fitting 


Experiment 


6s — 6p 
7 s — Qp 
7s — 7p 


6.464 
5.405 
13.483 


6.093 5.499 
5.450 5.198 
13.376 12.602 


5.510 
5.165 
12.641 


5.512 
5.201 
12.612 


5.497(8) a 
5.185(27) b 
12.625(18) c 


a Ref. ||. 
b Ref. [29]. 
c Ref. |30j. 












TABLE V. Percental 
approximations . 


$e deviation from experiment of calculated radial integrals in different 


Transition 




S( 2 ) 

with fitting 


Percental 


$e deviation 
t 


t 

with fitting 


6s — 6p 
7 s — 6p 
7s — 7p 




0.04 
0.3 
-0.2 




0.2 
-0.4 
0.1 


0.3 
0.3 
-0.1 



TABLE VI. Calculations of the hyperfine structure of Cs in different approximations. In the 
last column the experimental values are listed. Units: MHz. 



State 


RHF 


TDHF 


S(2) 

with fitting 


t 


S 

with fitting 


Experiment 


6s 


1425.0 


1717.5 


2306.9 


2287.9 


2285.7 


2298.2 a 


7s 


391.6 


471.1 


544.4 


539.6 


540.2 


545.90(9) b 


6P1/2 


160.9 


200.3 


291.5 


296.4 


293.3 


291.89(8) c 


7P1/2 


57.6 


71.2 


94.3 


95.4 


94.6 


94.35 a 



a Ref. @. 
b Ref. |2|. 
c Ref. @. 
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TABLE VII. Percentage deviation from experiment of calculated hyperfine structure constants 



in different approximations. 



State 


S( 2 ) 

with fitting 


Percentage deviation 
E 


E 

with fitting 


6s 


0.4 


-0.4 


-0.5 


7s 


-0.3 


-1.2 


-1.0 


6pi/ 2 


-0.1 


1.5 


0.5 


7pi/ 2 


0.05 


1.1 


0.3 



TABLE VIII. Percentage deviation from experiment of calculated v /hfs(s)hfs(p) (we will denote 



this by s — p in the tables) in different approximations. 



v /hfs(s)hfs(p) 


S( 2 ) 

with fitting 


Percentage deviation 
E 


S 

with fitting 


6s — 6p 


0.1 


0.5 


-0.02 


6s - 7p 


0.2 


-0.3 


-0.2 


7s — 6p 


-0.2 


0.2 


-0.3 



TABLE IX. Values for Ep^c i n different approximations; units 10 n ieaB(— Qw/N). 





£( 2 ) with fitting 


S 


E with fitting 


Epnc 


0.898 


0.902 


0.901 
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FIGURES 
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Xr 

FIG. 1. Lowest-order diagrams for PNC. Solid line denotes the bound electron; cross is the 
weak interaction; and dashed line is the El field. 
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(a) (b) (c) 

FIG. 2. Examples of diagrams representing the polarization of the atomic core by external 
fields. The dashed loop is the Coulomb field. (The diagrams we have presented are lowest-order 
exchange diagrams; there are also direct diagrams.) In diagrams (a) and (b) the core is polarized 
by a single field. Diagram (c) corresponds to the polarization of the core by both fields. 
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FIG. 3. Second-order correlation diagrams. Dashed line is the Coulomb interaction. Loop is 
the polarization of the atomic core. 
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FIG. 4. Chaining of the self-energy operator. 



FIG. 5. Screening of the Coulomb interaction. 
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FIG. 6. Hole-particle interaction in the polarization operator. 
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FIG. 7. The electron self-energy operator with screening and hole-particle interaction included. 
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FIG. 8. Lowest-order correlation corrections to the PNC El transition amplitude. 
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FIG. 9. Small corrections to the PNC El transition amplitude: external field inside the correla- 
tion potential. In diagrams (a) the weak interaction is inside the correlation potential (e>£ denotes 
the change in S due to the weak interaction); this is known as the weak correlation potential. 
Diagrams (b) represent structural radiation (photon field inside the correlation potential). 



o i 

(a) (b) 

FIG. 10. High-frequency contribution to radiative corrections. Diagram (a) corresponds to the 
Uehling potential. Diagram (b) is the vertex correction. The double line represents the bound 
electron; the single solid line is the free electron; the Coulomb interaction is denoted by the dashed 
line; and the filled circle denotes the nucleus. 
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